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Abstract- The radiative energy emitted by large fires and the 
corresponding smoke aerosol loading are simultaneously 
measured from the MODIS sensor from both the Terra and 
Aqua satellites. Quantitative relationships between the rates of 
emission of fire radiative energy and smoke are being developed 
for different fire-prone regions of the globe. Preliminary results 
are presented. When fully developed, the system will enable the 
use of MODIS direct broadcast fire data for near real-time 
monitoring of fire strength and smoke emission as well as 
forecasting of fire progression and smoke dispersion, several 
hours to a few days in advance. 

I. Introduction 

Fires are formidable agents of disaster, voraciously 
consuming natural and man-made resources and concurrently 
emitting particulate and gaseous pollutants into the 
atmosphere. Remote sensing provides a powerful means of 
fire monitoring and disaster mitigation. Each of the two Earth 
Observing System (EOS) Satellites, Terra and Aqua, is 
carrying the Moderate Resolution Imaging Spectroradiometer 
(MODIS) sensor aboard. MODIS has 36 spectral c hann els 
spread across file 0.4 to 14.4 pm wavelength range 
approximately, with at-nadir spatial resolutions of 0.25 km (2 
channels), 0.50 km (5 channels), and 1km (29 channels). 
Each MODIS covers nearly the entire globe once in the 
daytime and once in the nighttime daily, makin g a (Terra + 
Aqua) total of four measurements every 24 hours. 

Earlier fire remote sensing activities were limited to 
identifying fire locations, because the sensors used were 
easily saturated by strong signal from large or intense fires. 
With the enhanced radiometric range of file MODIS 4-pm 
fire channel, it is now possible to measure from space the rate 
of emission of radiative energy from fires [1]. In addition, 
MODIS operationally provides the total col umn loading and 
size parameters of smoke emitted by fires [2], among other 
aerosols. The availability of high quality fire and aerosol 
parameters from the same sensor, in an operational mode, 
offers a great opportunity for simultaneous monitoring of the 
severity of fire disasters and associated smoke pollution. 

Therefore, the focus of this application is to take 
advantage of this new data resource in developing tools for 
fire disaster management and air quality monitoring. From 
various field and airborne experiments, the rate of emission 
of fire radiative energy has been shown to be related to the 


rate of biomass consumption [3] as well as to the rate of 
smoke emission [1], The objective of this research is to 
establish this relationship between the rates of emission of 
radiative energy and smoke from satellite observations for 
fires in different biomes and regions. The derived biomass 
emission coefficients, together with fire radiative energy from 
the MODIS data acquired by the Forest Service Direct 
Broadcast system, will be assimilated into fire progression 
and Weather Research and Forecasting (WRF)-Smoke 
Dispersion models. This project will provide near real-time 
fire intensity and smoke emissions, as well as forecasts of fire 
spread rate and smoke dispersion, for use in decision support 
systems for wide-ranging applications, including for 
firefighting, alerts, early warning systems, evacuation, 
emissions inventory, weather, and climate. 

E. MODIS ACTIVE FIRE DATA 

The operational fire detection with MODIS is 
accomplished at 1-km spatial resolution, using mainly the 4- 
pm and 11-pm channels [1,4], However, the 0.66- pm and 
0.87-pm channels are used for filtering out pixels 
contaminated by sun glint. On detecting the 1-km resolution 
pixels containing fire, the MODIS fire algorithm uses the 4- 
pm brightness temperature of the fire pixel and the 
background to compute the rate at which the fire radiative 
energy (FRE) is emitted. This rate of emission of FRE per 
second (MJ/s), which is referred to as the Fire Radiative 
Power (FRP) in the MODIS fire product, is expressed in MW 
per pixel. The full description of the current MODIS fire 
algorithm, products, and validation activity is given in [4], 
while the detailed formulation for deriving FRP is given in 
[!]• 

m. Modis Aerosol Data 

Atmospheric aerosols include smoke, dust, sea salt, 
urban/industrial pollution, and other species. Aerosol 
parameters are retrieved from MODIS data at 10-km spatial 
resolution (level 2) only during the daytime, both over land 
and ocean [2]. Over land, the MODIS-retrieved aerosol 
parameters include: total column aerosol optical thickness 
(AOT or Tax) at 0.47, 0.55 and 0.66 pm wavelengths, 
Angstrom exponent (a), and the fraction of AOT contributed 


by the fine size mode aerosols (under the assumption that 
aerosols are composed of two main size modes: coarse and 
fine). AOT expresses the total column aerosol loading in each 
pixel, and is the most significant parameter used for various 
kinds of studies involving atmospheric aerosols. Although the 
fine mode fraction gives a general idea of aerosol type, 
MODIS does not specifically distinguish one aerosol type 
from another. Therefore, this study combining fire and 
aerosol data will enable accurate quantification of smoke 
aerosol, as distinct from other species. 

IV. MODIS FIRE/AEROSOL DATA ANALYSIS 

The basis of the application proposed in this paper is the 
relationship between the rates of emission of fire radiative 
energy and the corresponding smoke aerosol. To establish 
this relationship in a rudimentary way, various large regions 
of the world where fires occur seasonally were identified. 
Fig. 1 shows six regions on which the study was conducted, 
namely Southern Africa, Siberia, Western US, the Amazon, 
Northern Australia, and Western Russia. For each of the six 
regions, MODIS-measured FRP and AOT parameters were 
extracted for the year 2002 both from Terra- and Aqua- 
MODIS. 
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Fig. 1. Global map showing six selected regions used for the 
preliminary study of the relationships and cycles of rate of emission 
of fire radiative energy, Ej ■, and aerosol optical thickness, AOT, both 
derived from MODIS. 

Based on the assumption that the average lifetime of 
smoke in the atmosphere is about one week, weekly regional 
statistics of the FRP and AOT were computed for each of the 
regions for further analysis. For FRP, the weekly total was 
calculated for each region, since it is the total energy for the 
region (as opposed to average energy per pixel) that 
determines the emitted smoke amount However, to make the 
values comparable between regions, the total regional FRP is 
divided by the area of the region in deg 2 . On the other hand, 
since the smoke density would be expressed by the average 
AOT over each region, the weekly mean AOT was simply 
derived for each region. The total weekly FRP/deg 2 


(designated by Ej) was plotted against the mean weekly AOT 
for each region (Fig. 2a,b). There is an appreciably linear 
relationship between iyand the weekly mean AOT. Note that 
in i i.u me abscissa (the Ej axis) Fas a logarithmic scale. 


Therefore, for the three regions represented on Fig. 2b 
(Amazon, Australia, WRussia) AOT appears to be linearly 
related (not directly to Ej) but to log Ef. The relationships as 
well as other statistical information are listed in Table 1. It is 
important to emphasize that these results are only 
preliminary. Nevertheless these relationships have provided 
the encouragement to undertake more in-depth investigation 
in this area of research. 



Fig. 5. MODIS regional weekly mean aerosol optical thickness 
(AOT) at 550 nm against weekly total fire radiative power ( Ej) per 
degree squared (in GW). Best linear fits are shown for the group of 
points representing each region. Note that the horizontal axis on the 
left panel is plotted with a log scale. The correlation coefficients are 
listed in Table 1. 

V. RESEARCH DIRECTION AND APPLICATIONS 

Currently, research is in progress to develop clear and 
accurate relationships between the rates of emission of fire 
radiative energy and smoke aerosol. Traditional methods for 
deriving emission rates employ the use of emission factors e x 
(in g of emitted species x per kg of biomass burned), which is 
then used to estimate emission based on the mass of 
combusted biomass. However, the mass of burned biomass is 



difficult to measure from satellites. Current remote sensing 
methods attempt to estimate burned biomass mass either from 
the burned area (which is available only after the fire is over) 
or from fire pixel counts (which do not have a physical 
meaning, and can be confusing when the pixel spatial 
resolution changes). As a result, it has hitherto been difficult 
to monitor fire progression and smoke emission effectively in 
near real time. 

In this research, a new technique, using the fire radiative 
power (FRP) and aerosol optical thickness (both of which are 
measured by MODIS onboard the EOS Terra and Aqua 
satellites) to calculate FRE-based emission rates, is being 
developed. Therefore, instead of the traditional emission 
factors e x , a new parameter called “fire radiative energy 
emission coefficient”, C e (in g of emitted matter per Mjoule 
of radiative energy), which is related to the rate of biomass 
consumption, will be defined and derived for different 
regions and biomes. C e will be used in estimating smoke 
emissions from remotely-sensed FRP, which has the same 
physical meaning regardless of sensor or scale of observation. 

There is tremendous application potential for the 
availability of C e in combination with FRP, especially when 
received through the MODIS direct broadcast system in real 
time. FRP can be used to monitor the strength of fire in near- 
real time and can be multiplied by C e to derive the rate of 
generation of smoke, also in near-real time. Furthermore, this 


parameter (C e ) and measurement (FRP) can be applied in 
forecast models to predict the progression of fire and smoke 
dispersion for several hours or days in advance. Such 
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evacuation operations as well as for smoke pollution advance 
warning, in order to take measures to limit health risks. The 
results will also be very useful for keeping accurate inventory 
of fires and generated smoke, and for quantifying their effects 
on the environment, weather, and climate. 
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Table 1. Summary of the regional fire activity and smoke aerosol optical thickness for six biomass burning regions: For each 
region the table gives (from left to right): the region identification; Nf- the number of fire pixels detected by MODIS on 
Terra and Aqua during (July - December 2002); Ef ! pixel - the average rate of release of radiative energy per fire pixel 
in Mwatt; E/ deg 2 - the average rate of release of radiative energy per unit area of 1° x 1° latitude and longitude; AOT 
(r„) - the average regional aerosol optical thickness. The last two columns give the linear relationship and correlation, 
r, in each of the regions between the weekly total fire radiative energy and weekly average aerosol optical thickness 
throughout the fire season. Note that despite the variation of meteorological conditions that can deteriorate the 
relationship of AOT to fire activity, the correlation is still high -0.8. It is expected that once the wind field variation is 
accounted for in this research, the correlation will be higher. 


Regions 

N f 

£ ,/pixel (MW) 

Total £,/deg 2 

AOT (t a ) 

Best fit equation 

r 

S.Africa 

140413 

30 

27343 

0.26 

r,=0187+0.408F„ 

0.83 

W.Russia 

24858 

34 

2473 

0.31 

r a =0.37+0.078log(F w ) 

0.82 

Amazon 

137994 

40 

39527 

0.26 

r^O^+O.IOIogfF*,) 

0.67 

N. Australia 

86868 

48 

58807 

0.23 

t , =0.26+0.092log(F w ) 

0.93 

Siberia 

25938 

51 

3484 

0.25 

r a =0148+0.678F„ 

0.69 

US_West 

7720 

67 

529 

0.20 

r a =0141+0.888F w 

0.78 


N,= Number of fire pixels; £ f rate of emission of fire radiative energy; F w = weekly total £ f per degree squared 





